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Abstract. Introduced in the late 1960's [231113], the asymmetric exclusion pro- 
cess (ASEP) is an important model from statistical mechanics which describes a 
system of interacting particles hopping left and right on a one-dimensional lattice 
of n sites. It has been cited as a model for traffic flow and protein synthesis. In 
the most general form of the ASEP with open boundaries, particles may enter and 
exit at the left with probabilities a and 7, and they may exit and enter at the right 
with probabilities /3 and 5. In the bulk, the probability of hopping left is q times 
the probability of hopping right. The first main result of this paper is a combi- 
natorial formula for the stationary distribution of the ASEP with all parameters 
general, in terms of a new class of tableaux which we call staircase tableaux. This 
generalizes our previous work [141 115j for the ASEP with parameters 7 = 5 = 0. 
Using our first result and also results of Uchiyama-Sasamoto-Wadati [48] , we derive 
our second main result: a combinatorial formula for the moments of Askey- Wilson 
polynomials. Since the early 1980's there has been a great deal of work giving 
combinatorial formulas for moments of various other classical orthogonal polyno- 
mials (e.g. Hermite, Charlier, Laguerre, Meixner). However, this is the first such 
formula for the Askey- Wilson polynomials, which are at the top of the hierarchy 
of classical orthogonal polynomials. 



Contents 

1. Introduction 

2. The asymmetric exclusion process (ASEP) 

3. Staircase tableaux and the stationary distribution of the ASEP 

4. Askey- Wilson polynomials and a formula for their moments 

5. A more flexible Matrix Ansatz 

6. The proof of the stationary distribution 

7. The proof of our Askey- Wilson moment formula 

8. Open problems 

9. Appendix: Staircase, permutation, and alternative tableaux 
References 



10 
I I 

20 
21 

22 
23 



2000 Mathematics Subject Classification. Primary 05E10; Secondary 82B23, 60C05. 

Key words and phrases, permutation tableaux, asymmetric exclusion process, Matrix Ansatz, 
staircase tableaux, Askey- Wilson polynomial. 

Both authors were partially supported by the grants ANR blanc Gamma and ANR-08-JCJC- 
0011, and the second author was partially supported by the NSF grant DMS-0854432 and an Alfred 
Sloan Fellowship. 

1 



2 



SYLVIE CORTEEL AND LAUREN K. WILLIAMS 



1. Introduction 

The asymmetric exclusion process (ASEP) is an important model from statistical 
mechanics which was introduced independently in the context of biology [33] and in 
mathematics [13] around 1970. Since then there has been a huge amount of activity 
on the ASEP and its variants for a number of reasons: although the definition of 
the model is quite simple, the ASEP is surprisingly rich. For example, it exhibits 
boundary-induced phase transitions, spontaneous symmetry breaking, and phase 
separation. Furthermore, the ASEP is regarded as a primitive model for translation 
in protein synthesis [33J, traffic flow [39], and formation of shocks [19]; it also appears 
in a kind of sequence alignment problem in computational biology [91]. Important 
mathematical techniques used to understand this model include the Matrix Ansatz 
of Derrida, Evans, Hakim, and Pasquier [18] . and also the Bethe Ansatz [16l |26| 142] . 

This paper concerns the ASEP on a one-dimensional lattice of n sites with open 
boundaries. Particles may enter from the left at rate adt and from the right at rate 
5dt; they may exit the system to the right at rate (3dt and to the left at rate ^dt. 
The probability of hopping left and right is qdt and udt, respectively^ 

Dating back to at least 1982, it was realized that there were connections between 
this model and combinatorics - for example, Young diagrams appeared in [30], and 
Catalan numbers arose in the analyses of the stationary distribution of the ASEP in 
P21 HE]. More recently, beginning with papers of Brak and Essam [7], and Duchi and 
Schaeffer [20] , and followed by a number of works including pfl [El EED, El EE EQ] , there 
has been a push to understand the roots of the combinatorial features of the ASEP. 
Ideally, the goal is to find a combinatorial description of the stationary distribution: 
that is, to express each component of the stationary distribution as a generating 
function for a set of combinatorial objects. Up to now, the best result in this direction 
was provided by [HI [15], which addressed this question when 7 = 5 = in terms 
of permutation tableaux, combinatorial objects introduced in the context of total 
positivity on the Grassmannian [36] . 

There is a second reason why combinatorialists became intrigued by the ASEP. 
Papers of Sasamoto [38] and subsequently Uchiyama, Sasamoto, and Wadati [18] 
linked the ASEP with open boundaries to orthogonal polynomials, in particular, 
to the Askey- Wilson polynomials. The Askey- Wilson polynomials are orthogonal 
polynomials that sit at the top of the hierarchy of (basic) hypergeometric orthogonal 
polynomials [2j |25j [32], in the sense that all other polynomials in this hierarchy 
are limiting cases or specializations of the Askey- Wilson polynomials. It is known 
that orthogonal polynomials have many combinatorial features: indeed, starting 
around the early 1980's, mathematicians including Flajolet [23], Viennot [19], and 
Foata [21], initiated a combinatorial approach to orthogonal polynomials. Since then, 
combinatorial formulas have been given for the moments of (the weight functions of) 



Actually there is no loss of generality in setting u — 1, so we will often do so. 
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many of the polynomials in the Askey scheme, including g-Hermite, Chebyshev, q- 
Laguerre, Charlier, Meixner, and Al-Salam-Chihara polynomials, see e.g. [271 I2H1 
130] |2U HI]. However, no such formula was known for the moments of the Askey- 
Wilson polynomials. Therefore when the paper [48J provided a close link between 
the moments of the Askey- Wilson polynomials and the stationary distribution of the 
ASEP, this gave the idea that a complete combinatorial understanding of the ASEP 
might also give rise to the combinatorics of the Askey- Wilson moments. 

In this paper we give a complete solution to both of the above problems. Namely, 
we introduce some new combinatorial objects, which we call staircase tableaux, and 
we prove that generating functions for staircase tableaux describe the stationary 
distribution of the ASEP, with all parameters general. We then use this result 
together with jH] to give a combinatorial formula for the moments of the Askey- 
Wilson polynomials. 

The method of proof for our stationary distribution result builds upon impor- 
tant work of Derrida, Evans, Hakim, and Pasquier [18J, who introduced a Matrix 
Ansatz as a tool for understanding its stationary distribution. Briefly, the Matrix 
Ansatz says that if one can find matrices and vectors satisfying certain relations (the 
DEEP algebra), then each component of the stationary distribution of the ASEP can 
be expressed in terms of certain products of these matrices and vectors. Knowing 
this Ansatz, the strategy for proving our stationary distribution result which one 
would like to employ is the following: find matrices and vectors satisfying the DEHP 
algebra, and show that appropriate products enumerate staircase tableaux. 

However, it has been known since 1996 [211 Section IV] that when the parameters 
of the ASEP satisfy certain relations (for example afi = q % ^5), there is no represen- 
tation of the DEHP algebra. Therefore the above strategy cannot succeed when all 
parameters of the ASEP are general. What we do instead is to introduce a slight 
generalization of the Matrix Ansatz (see Theorem 15. 2p . which is more flexible, albeit 
harder to use: instead of checking three identities, one must check three infinite fami- 
lies of identites. So our strategy is to find matrices and vectors such that appropriate 
products enumerate staircase tableaux, and prove that they satisfy the relations of 
the Generalized Matrix Ansatz. This second step is quite involved, as our "matrices" 
and "vectors" are somewhat complicated (they have four and two indices each, see 
Section I6.1D . and there is no easy way to use induction to prove the three infinite 
families of identities. 

We believe that our new staircase tableaux deserve further study, because of their 
combinatorial interest and their potential connection to geometry. For example, 
staircase tableaux of size n have cardinality A n n\, and hence are in bijection with 
doubly-signed permutations. In [13], we will prove this with an explicit bijection, 
and also develop connections to other combinatorial objects, including matchings, 
permutations, and trees. Furthermore, because of the connection to the ASEP, we 
know that our staircase tableaux have some hidden symmetries which are not at 
all apparent from their definition. For instance, it is clear from the definition of 
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the ASEP that the model remains unchanged if we reflect it over the y-axis, and 
exchange parameters a and 5, (3 and 7, and q and u. However, the corresponding 
bijection on the level of tableaux has so far eluded us. Finally, staircase tableaux 
generalize permutation tableaux, which index certain cells in the non-negative part of 
the Grassmannian [36J; it would be interesting to better understand the relationship 
between the tableaux, the ASEP, and the geometry, and potentially generalize it to 
staircase tableaux. 

It is worth mentioning that in recent years there has been an explosion of activity 
[3, El E21 El ETl SI S5J SSI ST] surrounding another version of the ASEP, in which 
particles hop not on a finite lattice but on Z. Much of this interest was inspired 
by Johansson's work [29], which reinterprets the TASEP (totally asymmetric ex- 
clusion process) as a randomly growing Young diagram. Johansson then used the 
well-understood combinatorics of Young diagrams and semi-standard tableaux to 
approach the problem of current fluctuations in the TASEP. In light of this, one may 
hope that a better understanding of the combinatorics of staircase tableaux could 
lead to even more results on the ASEP with open boundaries. 

The structure of this paper is as follows. In Section [2] we define the ASEP (with 
open boundaries), and in Sections [3] and 0] we state our main results on the ASEP and 
on Askey- Wilson polynomials. In Section [5] we prove a generalized Matrix Ansatz, 
and in Sections [6] and [7] we prove our results on the ASEP and Askey- Wilson poly- 
nomials. Section [8] gives open problems, and the Appendix describes how staircase 
tableaux generalize permutation tableaux and alternative tableaux [50J. 

Acknowledgments: We would like to thank I. Gessel, R. Stanley, and J. de 
Gier for interesting remarks. We are also grateful to the referees for their thoughtful 
comments, which helped us to greatly improve the exposition. 

2. The asymmetric exclusion process (ASEP) 

The ASEP is often defined using a continuous time parameter [18]. However, one 
can also define it as a discrete-time Markov chain, as we will do below. The contin- 
uous and discrete-time definitions are equivalent in the sense that their stationary 
distributions are the same |20j. 

Definition 2.1. Let a, (3, 7, 5, q, andu be constants such that 0<a<l,0</3<l ; 
< 7 < 1, < 5 < 1, < q < 1, and < u < 1. Let B n be the set of all 2 n 
words in the language {o, •}*. The ASEP is the Markov chain on B n with transition 
probabilities: 

• If X = A • oB and Y = A o »B then Px,y = (particle hops right) and 
Py,x — (particle hops left). 

• If X = oB and Y = »B then Pxy — ^xj (particle enters from the left). 

• If X = B» and Y = Bo then Px,y — ^rj (particle exits to the right). 

• If X = »B and Y = oB then Px,y = ^fj (particle exits to the left). 

• If X = Bo and Y = B» then Px,y = ^tt (particle enters from the right). 
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• Otherwise Px,y = for Y ^ X and Px,x = 1 — TIx^y Px,y- 

Note that we will sometimes denote a state of the ASEP as a word in {0, l} n and 
sometimes as a word in {o, • }". In these notations, the symbols 1 and • denote a 
particle, while and o denote the absence of a particle, which one can also think of 
white particle. 

See Figure [1] for an illustration of the four states, with transition probabilities, for 
the case n = 2. The probabilities on the loops are determined by the fact that the 
sum of the probabilities on all outgoing arrows from a given state must be 1. 




Figure 1 . The state diagram of the ASEP for n = 2 

In the long time limit, the system reaches a steady state where all the probabilities 
Pn(Ti, Tii ■ ■ ■ i T n) of finding the system in configurations (n, r 2 , . . . , r n ) are stationary. 
More specifically, the stationary distribution is the unique (up to scaling) eigenvector 
of the transition matrix of the Markov chain with eigenvalue 1. 

The ASEP clearly has multiple symmetries, including the following. 

• The "left-right" symmetry: if we reflect the ASEP over the y-axis, we get 
back the same model, except that the parameters a and 5, 7 and 0, and u 
and q are switched. 

• The "arrow-reversal" symmetry: if we exchange black and white particles, 
we get back the same model, except that the parameters a and 7, (3 and S, 
and u and q are switched. 

• The "particle-hole" symmetry: if we compose the above two symmetries, i.e. 
reflect the ASEP over the y-axis and exchange black and white particles, we 
get back the same model, except that a and (3, and 7 and 5 are switched. 

These symmetries imply results about the stationary distribution. 

Observation 2.2. The steady state probabilities satisfy the following identities: 

• Pn(T~i, . . . , Tn) = P n (T n , . . . , T\ ) | a ^S t p^^ )U ^q, 

• Pn\Tli ■ ■ ■ 1 Tn) = PniX 7*1) ■ ■ ■ j 1 Tn) \ 0-0-7, /9<-^(5,w«-<?; 

• Pn{Tl, ■ ■ ■ , T n ) = P n (l — T n , . . . ,1 — Tl)\ a ^fi^S- 

Above, the notation \ a ++g indicates that the parameters a and 5 are exchanged. 
These symmetries are related to the symmetries of the Askey- Wilson polynomials 
(see for example Remark 14.1 p . though neither is a direct consequence of the other. 
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3. Staircase tableaux and the stationary distribution of the ASEP 

The main combinatorial objects of this paper are some new tableaux which we call 
staircase tableaux. These tableaux generalize permutation tableaux (equivalently, 
alternative tableaux). 

Definition 3.1. A staircase tableau of size n is a Young diagram of "staircase" 
shape (n, n — 1, . . . , 2, 1) such that boxes are either empty or labeled with a, (3, 7, or 
5, subject to the following conditions: 

• no box along the diagonal is empty; 

• all boxes in the same row and to the left of a (3 or a 5 are empty; 

• all boxes in the same column and above an a or a 7 are empty. 

The type of a staircase tableau is a word in {o, •}" obtained by reading the diagonal 
boxes from northeast to southwest and writing a • for each a or 5, and a o for each 
(3 or 7. 

See the left of Figure [2] for an example of a staircase tableau. 
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Figure 2. A staircase tableau of size 7 and type o o • • • o o 

Definition 3.2. The weight wt(T) of a staircase tableau T is a monomial in a, (3, 7, S, q, 

and u, which we obtain as follows. Every blank box of T is assigned a q or u, based 
on the label of the closest labeled box to its right in the same row and the label of the 
closest labeled box below it in the same column, such that: 

• every blank box which sees a (3 to its right gets assigned a u; 

• every blank box which sees a 5 to its right gets assigned a q; 

• every blank box which sees an a or 7 to its right, and an a or 5 below it, gets 
assigned a u; 

• every blank box which sees an a or 7 to its right, and a /3 or 7 below it, gets 
assigned a q. 

After assigning aqoruto each blank box in this way, the weight ofT is then defined 
as the product of all labels in all boxes. 

The right of Figure [2] shows that this staircase tableau has weig ht a 3 /3 2 

Remark 3.3. The weight of a staircase tableau always has degree n(n + l)/2. For 

convenience, we will sometimes set u = 1, since this results in no loss of information. 
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Our first main result (to be proved in Section [6]) is the following. 

Theorem 3.4. Consider any state r of the ASEP with n sites, where the parameters 
a, (3, 7, 5, q and u are general. Set Z n = J2t w ^0~)> where the sum is over all staircase 
tableaux of size n. Then the steady state probability that the ASEP is at state r is 
precisely 

E r wt(T) 



where the sum is over all staircase tableaux T of type r. In particular, Z n is the 
partition function for the ASEP. 



Figure [3] illustrates Theorem 13. 41 for the state •• of the ASEP. All staircase tableaux 
T of type •• are shown. It follows that the steady state probability of •• is 

a 2 u + S 2 q + a5q + a5u + a 2 5 + a(35 + a^5 + aS 2 
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Figure 3. The tableaux of type •• 

We can also obtain combinatorial formulas for various physical quantities. Theo- 
rem 13.51 below will be proved in Section 16.41 

Theorem 3.5. Consider the ASEP with n sites. Then we have the following: 

• The current in the steady state is Zn ~ l( -" / ^~ 7<59 — -, where Z n is the generating 
function for staircase tableaux of size n. 

• The average particle number (rj) n at site i is given by Z~ x times the generating 
function for all staircase tableaux of size n which have an a or 5 at the ith 
position along the diagonal. 

• Similarly, the m-point function (r^ . . -T im ) n is given by Z~ x times the gen- 
erating function for all staircase tableaux of size n which have an a or 5 in 
positions i 1 ,i 2 , . . . , and i m along the diagonal. 

Remark 3.6. In [H] and [15] , which concerned the ASEP with parameters 7 = 
5 = 0, we gave combinatorial expressions for the stationary distribution in terms of 
permutation tableaux: more specifically, the steady state probability of a state r was 
given by enumerating permutation tableaux lying in a Young diagram of shape A(r) ; 
according to the number of 1 's in the top row and the number of unrestricted rows. 

Subsequently [10] and [12] observed that a permutation tableau is determined by 
complementary statistics, namely the positions of the topmost 1 's not in the first 
row, and rightmost restricted zero's. Independently, Viennot [50] defined some new 
alternative tableaux, in order to give a simpler formula for the stationary distribution 
of the ASEP when 7 = 5 = than the one in [HI [15j ; it turns out that alternative 
tableaux can be obtained from permutation tableaux, by making boxes red, blue, or 
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empty, based on whether the box of the associated permutation tableau contains a 
topmost 1, rightmost restricted zero, or neither. Our staircase tableaux generalize 
both kinds of tableaux, but now have non-empty boxes labeled by a,/3,j,5. Addi- 
tionally, instead of working with Young diagrams of various shapes, we now always 
work with a staircase shape, whose type encodes the same information (namely, the 
corresponding state of the ASEP) that the shape used to. Working with this larger 
shape seems to be the only natural way to assign the appropriate powers of q and u 
to the tableau. 



Remark 3.7. .4s an alternative to Definition ^. £X suppose we define the dual weight 
wt'(T) of a staircase tableau by taking the product of all labels in all boxes, after 
having filled the blank boxes of T according to the following rule: 

• every blank box which sees an a below it gets assigned a u; 

• every blank box which sees a 7 below it gets assigned a q; 

• every blank box which sees an a or 5 to its right, and a (3 or 5 below it, gets 
assigned a q; 

• every blank box which sees a (3 or 7 to its right, and a (5 or 5 below it, gets 
assigned a u. 



Then Theorem \ 3.4\ continues to hold, with wt replaced by wt'. This follows from 
the left-right symmetry of the ASEP. More specifically, note that if T is a staircase 
tableau, then the tableau T' obtained by transposing T then switching a and 5, and 



(3 and 7 is still a staircase tableau. Our observation now follows from Theorem\3.4 



the fact that wt'(T') = wt(T), and Observation \2.S[ (Alternatively, we could have 



proved the observation using a method analogous to the proof of Theorem \3.4\ ) 



4. Askey- Wilson polynomials and a formula for their moments 

The Askey- Wilson polynomials are orthogonal polynomials with five free parame- 
ters (a, b, c, d, q). They reside at the top of the hierarchy of the one- variable orthog- 
onal polynomial family in the Askey scheme [2J [25j [32] . In this section we define 
the Askey- Wilson polynomials, following the exposition of [2] and [48] , then state a 
combinatorial formula for their moments. 

The g-shifted factorial is defined by 

s n— 1 

(a u a 2 , ■ ■ • , a s ; q) n = ] [ - a r q k ), 

r=l k=0 

and the basic hypergeometric function is given by 



at, 



,q,z 



E 

fc=0 



[ai, ■ ■ ■ ,a r ;q)k ^_^k^k(k-i)/2y+s-r z k 



(&!,••• ,b s ,q;q) k 



The Askey- Wilson polynomial P n (x) 
P n ( x ) — a~ n (ab, ac, ad; q) n 4 3 



= P n (x; a, b, c, d\q) is explicitly defined by 
q~ n , q n ~ 1 abcd, ae t9 , ae~ t6 



ab, ac, ad 
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with x = cos8 for n G Z + := {0, 1, 2, ■ ■ ■ }. It satisfies the three-term recurrence 

A n P n+1 (x) + B n P n (x) + C n P n „i(x) = 2xP n (x), 
with Pq(x) = 1 and P_i(x) = 0, where 
I — q n ~ 1 abcd 



A 



[l - q 2n - l abcd)(l - q 2n abcd) ' 



n-l 

_ (I — g n )(l - g n - 1 a6)(l - q n - l ac)(l - q n - l ad)(l - ^^(l - g n - x 6d)(l - q n ~ x cd) 
n ~ (f - q 2n - 2 abcd)(l - q 2n - l abcd) ' 

and s = a + b + c + d, s' = a+ x + b^ 1 + c+ x + d~ . 

Remark 4.1. It is obvious from the three-term recurrence that the polynomials P n {x) 
are symmetric in a, b, c and d. 

For \a\, \b\, \c\, \d\ < 1, using z = e %e , the orthogonality is expressed by 

dz ( z + z- l \ ( z + z- l \ ( z + z~ x \ 

w I - I P m P n = h n d mn , 



c Aniz 



where the integral contour C is a closed path which encloses the poles at z — aq k , bq k , 
cq k , dq k (k G Z+) and excludes the poles at z — (aq k )~ l , (cg fc ) _1 , (<ig fe )~ 1 

(k G Z + ), and where 



w(cosB) - 



(ae ie , ae~ ie , be ie , 6e" ie , ce ie , ce~ ie , de ie , de~ ie ; q) 
h n (I — q n ~ 1 abcd)(q ) ab, ac, ad, be, bd, cd; q) n 
ho (1 — q 2 ™" 1 abed) (abed; q) n 

(abed; q)^ 



oo 



/?() 



(g, ab, ac, ad, be, bd, cd; q)^ 



(In the other parameter region, the orthogonality is continued analytically.} 
The moments are defined by 

dz ( z + z~ x \ ( z + z~ } x '' 



l c 4mz 



The second main result of this paper is a combinatorial formula for the moments 
of the Askey- Wilson polynomials. In Theorem 14.21 below, we use the substitution 

1 - Q a 1 ~ a 

a = , p 



1 + ac + a + c 1 + bd + b + d ' 

— (1 — q)ac — (l — q)bd 

7 = i~, ; — , o 



1 + ac + a + c 1 + bd + b + d 
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which can be inverted via 



b = 
d = 
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Recall that Zi = ^t w ^(^~)> wriere the sum is over all staircase tableaux of size I. 
Theorem 4.2. The kth moment of the Askey- Wilson polynomials is given by 

5. A MORE FLEXIBLE MATRIX ANSATZ 

One of the most powerful techniques for studying the ASEP is the so-called Matrix 
Ansatz, an Ansatz given by Derrida, Evans, Hakim, and Pasquier [18] as a tool for 
solving for the steady state probabilities P u (ti, . . . , r n ) of the ASEP. In this section 
we will start by recalling their Matrix Ansatz, and then give a new generalization of 
it which we require for our proof of Theorem 13.41 

For convenience, in this section we set u = 1. Also, we define unnormalized weights 
fniji, . . . , r n ), which are equal to the P n {Ti, . . . , r n ) up to a constant: 

where Z n is the partition function fn{ T i, ■ ■ ■ > T n)- The sum defining Z n is over all 
possible configurations r G {0, l} ra . Derrida et al showed the following. 

Theorem 5.1. [18J Suppose that D and E are matrices, V is a column vector, and 
W is a row vector, with WV = 1, such that the following conditions hold: 

DE - qED = D + E, (3DV - 5EV = V, aWE - yWD = W. 

Then for any state r = (ti, . . . , r n ) of the ASEP, 

n 

f n (n, . . . , O = wQJfcD + (l - n)E))V. 

Note that nr=i( r *-^ + 0- ~~ T i)E) is simply a product of n matrices D or E with 
matrix D at position i if site i is occupied (t{ = 1). Also note that Theorem 15.11 
implies that Z n = W(D + E) n V . 

We now state and prove a more flexible version of Theorem 15.11 Our proof gener- 
alizes the argument given in [18] . 
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Theorem 5.2. Let {A n } n >o be a family of constants. Let W and V be row and 

column vectors, with WV = 1, and let D and E be matrices such that for any words 
X and Y in D and E, we have: 

(I) WX(DE - qED)YV = X ix \+\y\+2W X (D + E)YV ; 
(II) (5WXDV - 5WXEV = \\ X \+iWXV; 
(III) aWEYV - jWDYV = X lY]+1 WYV. 

(Here \X\ is the length of X .) Then for any state r = (n, . . . , r n ) of the ASEP, 

n 

(1) /n(r) = W(jJ(riD + (1 - n)E))V. 

i=l 

Proof. We are in the steady state of the ASEP if the net rate of entering each state 
(ti, . . . , t„) is 0, or in other words, the following expression equals 0: 

(2) (-ir (-af n (0, r 2 , . . . , r n ) + 7 /n(l, r 2 , . . . , r n )) 

n-l 

( 3 ) + ^(-Wxfa ^ Ti+MUfa, . . . , 1, o, . . . , O - g/ n (n, . . . , 0, 1, . . . , r n )) 

i=l 

(4) +(_l)^(/3/ n ( ri , . . . , r n _ 1; 1) - 5/ n ( Tl , . . . , r n _ 1; 0)). 

In ([3]) above, the arguments 1, and 0, 1 are in positions i and i + 1, and x is the 
boolean function taking value 1 or based on whether its argument is true or false. 
So what we have to prove is that the quantities in the right-hand-side of equation 
([[]) satisfy this equation. 

By the assumptions (I), (II), and (III) of the theorem, we have the following: 

• the expression (j2J) equals ±A n / n _i(r 2 , . . . , r n ), based on whether n is 1 or 0; 

• (G2) equals when n = Tj+i; otherwise, based on whether n is 1 or 0, it equals 
T^E"=i(/"-i( T i' . . . ,Ti, . . . ,r„) + fn-iin, . . . ,fi+i, . . . ,r n )); 

• and (jlj) equals =FA n / n _i(ri, . . . , r n _i), based on whether r n is 1 or 0. 

(Here f% denotes the omission of the ith component.) Then using these conditions, 
it is easy to verify that the sum of (TSJ), (jHJ), and (jlj) is equal to 0, since all terms 
involving / n _x cancel out. 

□ 

6. The proof of the stationary distribution 

In this section we will prove Theorem 13.41 by: defining vectors W, V and matrices 
D, E; proving that they have the requisite combinatorial interpretation in terms of 
staircase tableaux; and checking that they satisfy the relations of Theorem 15. 2\ with 
Ao = 1 and A n = a/3 — 7<5g™~ 1 for n > 1. 

This is analogous to the proof of [TU Theorem 3.1], albeit much more difficult: in 
[T4] . it was obvious that our matrices and vectors satisfied the Matrix Ansatz [T4l 
Lemma 2.5], and easy to show that our combinatorial objects were described by the 
algebraic relations of the Ansatz, see [TU Figure 6] and the surrounding discussion. 
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In contrast, in this more general situation, we can give a combinatorial proof of 
relation (III) of our new Matrix Ansatz, but not for (I) or (II). Instead we give a 
rather difficult algebraic proof of (I) and (II). First of all, our new "vectors" and 
"matrices" have two and four indices, respectively, which makes working with them 
more complicated. Second, to use Theorem 15.21 instead of proving that our vectors 
and matrices satisfy three identities (as in Theorem 15. ip . we must prove that they 
satisfy three infinite families of identities. Moreover, there is no obvious way to use 
induction to prove these identities: one cannot take one of the identities and multiply 
on the left or right to obtain the next identity in the family. 

Remark 6.1. In this section we assume u — 1. Recall that this is no loss of gen- 
erality, as the weight of a staircase tableau of size n is always a monomial of degree 
n(n + l)/2. However, setting u = 1 has the advantage of making our bookkeeping a 
little easier. 

6.1. The definition of our matrices. 

Definition 6.2. In what follows, indices range over the non-negative integers. In 
particular, our matrices and vectors are not finite. We define "row" and "column" 
vectors W = (Wik)i,k and V = {Vje)j,e., and matrices D = (Dij^,e)i,j,k,e and E = 
,j,k,i)i,j,k,e by the following: 



W, 



k 



1 if i = k = 0, 
otherwise, 



Vj£ = 1 always. 



D 



i,j,k,l 



E 



>i,j,h,£ 



if j < i or £ > k + 1, 

5q l if i = j — 1 and k = £ = 0, 

aq % if i = j, k = and i = 1, 

5(Aj-i,fc-i,£ + #ij_i,jfe-i I *) + D idi k-i,e-i otherwise. 

if j < % or I > k + 1, 

(3q l if i = j and k = £ = 0, 

7g* if i — j, k — and £ — 1, 

(3(Di,j,k-i,i + E ijj)k _ 1>e ) + qEij >k _x^ x otherwise. 
By convention, if any subscript i, j, k or £ is negative, then D^ki — ^ijU = 0- 
Example 6.3. 

-^0,2, 2,0 = /3(-Do,2,l,0 + -^0,2,1,0) + 9-^0,2,1,-1 

= /3(-Do,2,i,o + -^0,2, 1,0) 

= (3[5(D 0t i ;0 ,o + -^o, 1,0,0) + Ai,2,o,-i + /?(A),2,o,o + -^0,2,0,0) + ^0,2,0,-1] 
= /3<5(-Do,i,o,o + -Eo,i,o,o) + /^ 2 (-Do,2,o,o + -^0,2,0,0) 
= 05(8 + 0) + /3 2 (0 + 0) = (35 2 
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Here, we think of the two coordinates i and k as specifying a "row" of a matrix, 
and the two coordinates j and £ as specifying a "column" of a matrix. Therefore 
matrix multiplication is defined by 

(MN) i!j>k>£ = J2 M ^,k,b N a,j,b,e- 

0,6 

Note that when M and N are products of .D's and E's the sum on the right-hand-side 
is finite. Specifically, if M is a word of length r in D and E, then & & is unless 
a + 6<z + A;-|-r. This can be shown by induction from the definition of D and E, 
or by using the combinatorial interpretation given in Lemma 16.51 

6.2. The combinatorial interpretation of our matrices in terms of tableaux. 

We say that a row of a staircase tableau T is indexed by (3 if the leftmost box in 
that row which is not occupied by a q or u is a (3. Note that every box to the left 
of that (3 must be a u. Similarly we will talk about rows which are indexed by 5; in 
this case, every box to the left of that 5 must be a q. We will also talk about rows 
which are indexed by a/7, which is shorthand for rows which are indexed by a or 7. 

Theorem 6.4. If X is a word in D 's and E 's, then: 

• Xijke is the generating function for all ways of adding \X\ new columns to a 
staircase tableau with i rows indexed by 5 and k rows indexed by a/7, so as 
to obtain a new tableau with j rows indexed by 5 and £ rows indexed by a/7. 

• (WX)ji is the generating function for staircase tableaux of type X which have 
j rows indexed by 5 and £ rows indexed by a/7 (and hence \X\ — j — £ rows 
indexed by (3.) 

• WXV is the generating function for all staircase tableaux of type X . 

The main step in proving Theorem 16.41 is the following lemma, which says that 
the matrices D and E are "transfer matrices" for building staircase tableaux. 

Lemma 6.5. Dij^,e is the generating function for the weights of all possible new 
columns with an a or 5 in the bottom box that we could add to the left of a staircase 
tableau with i columns indexed by S and k columns indexed by a or 7, obtaining a 
new staircase tableau which has j columns indexed by 5 and £ columns indexed by a 
or 7. Similarly for Eij^j, where the new column has a (3 or 7 in the bottom box. 

Proof. Let D[- M denote the generating function for all possible new columns with an 
a or 5 in the bottom box that we could add to the left of a staircase tableau with 
i columns indexed by 8 and k columns indexed by a/7, obtaining a new staircase 
tableau which has j columns indexed by 5 and £ columns indexed by a or 7. We will 
show that D[- M = E^i by showing that D' satisfies the same recurrences. 

Note that D[- M = if j < i because adding a new column to a staircase tableau 
never decreases the number of rows indexed by 5. Also D[- M = Oif^>A; + l because 
when we add a new column we can never increase the number of rows indexed by 
a/7 by more than 1. 
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Now suppose that k = 0. If we are starting from a tableau with i rows indexed 
by 5 and rows indexed by a/7, then the only way to add a new column is to add 
a column with an a or 5 at the bottom, with all boxes above empty. If we add a 5, 
then the resulting tableau has £ = rows indexed by a/7 and j — i + 1 rows indexed 
by 5. The weight of the new column will be 5q l . On the other hand, if we add an a 
at the bottom, then the resulting tableau has £ = 1 rows indexed by a/7 and j = i 
rows indexed by 5. The weight of the new column will be aq % . From this discussion 
it follows that D[- M = 5q l when j = i + 1 and k = £ = 0, and D[- M = aq 1 when 
j = i, k = 0, and £ = 1. 

In all other situations, we can assume that k > 1. Suppose that we are adding a 
new column C with an a or 5 at the bottom to the left of a staircase tableau with i 
columns indexed by S and k columns indexed by a/7, so as to create a new tableau 
T. Consider the lowest box B of C whose row in T is indexed by an a or 7 (such 
a box exists since k > 1). If we fill B with an a, 0, 7 or 5, then the bottom box of 
C must contain a 5. In this case, if we ignore that bottom S, then our choices for C 
are exactly the same as our choices would be for adding a new column to the left of 
a staircase tableau with i columns indexed by S and k — 1 columns indexed by a/7. 
Therefore, filling B with an a, /3, 7 or 5 gives us a contribution of d(D' + -E")m'-i.£i-i/ 
to our generating function. 

On the other hand, if we leave B empty, then this box will get a weight u = 1 
(recall Remark 16.11) . Filling the rest of the column C is like adding a new column to a 
staircase tableau with % rows indexed by 5 and k — 1 rows indexed by a/7. Therefore 
leaving B empty gives us a contribution of D' i j k _ i e _ 1 to our generating function. 

It follows that when k > 1, D' ijU = 5{D' + + D 'i,j,k-i,i-i- 

Similarly, we define E[- M to be the generating function for all possible new columns 
with a (3 or 7 in the bottom box that we could add to the left of a staircase tableau 
with i columns indexed by 8 and k columns indexed by a/7, obtaining a new staircase 
tableau which has j columns indexed by 5 and i columns indexed by a/7. The proof 
that E[- M = Eij k £ is analogous to the proof we gave for D' . □ 

Proof of Theorem \6.4\ The first item follows from Lemma [6.51 and the definition of 
matrix multiplication. Multiplying at the left by a W has the effect that we start 
with the empty tableau and then add columns according to X: so WXji is the 
generating function for staircase tableaux of type X which have j rows indexed by 
5 and £ rows indexed by a/7. Finally, multiplying WXj£ on the right by V has 
the effect of summing over all 5 and £, so WXV is the generating function for all 
staircase tableaux of type X. □ 

6.3. The proof that our matrices satisfy the Matrix Ansatz. We now prove 
that our matrices satisfy Theorem 15.21 with A n = a/3 — / -f5q n ~ 1 for n > 1. Relation 
(III) has a simple combinatorial proof. However, this proof does not work for relation 
(II), and indeed it will require a lot more work to prove (I) and (II). 

Lemma 6.6. Relation (III) of Theorem \5.!2\ holds. 
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Proof. Using Theorem I6.4[ relation (III) can be reformulated in terms of staircase 
tableaux. First we rewrite (III) as 

aWEYV + ~ibq n - x WYV = -fWDYV + afiWYV, 

where n—1 = \Y\. Since a "type E" corner box of a staircase tableau must be either 
a or 7, we can rewrite this again as 

aWEpYV + aWE^YV + -f5q n ' l WYV = ^WD a YV + jWD s YV + a/3WYV. 

Here WEpYV denotes the generating function for staircase tableaux of type EY, 
whose northeast corner box is a (3; the terms WE^YV, WD a YV, and WDsYV are 
defined analogously!^ 

It is now easy to see that aWEpYV = a(3WYV and ^WD a YV = ^5q n ~ l WYV , 
since a box labeled /3 must have only empty boxes (weighted u = 1) to its left, and 
a box labeled 6 must have only empty boxes (weighted q) to its left. Also, since 
the rules for the weight of an empty box which sees a 7 to its right are the same as 
the rules for the weight of an empty box which sees an a to its right, we have that 
aWE^YV = -fWD a YV. This proves relation (III). □ 

Lemma 6.7. For any word Y in D and E, we have Y^i = g' y ' Vj-i.j-i^- 

Proof. We use Theorem 16.41 Note that both Y^ and i^-ij-i^/ enumerate the ways 
of adding \Y\ new columns to a staircase tableau T so as to increase by j — i the 
number of rows indexed by 5, and to increase by I — k the number of rows indexed 
by a/7. The only difference is the initial number of rows indexed by 5. Since Y^g 
has one extra initial row indexed by 5, this will contribute \Y\ extra empty boxes 
which all get the weight q. Therefore Y^u = Q^Yi-ij-i^i- d 

Proposition 6.8. To prove (I) and (II), it suffices to prove the following identities 
for all non-negative integers j and I: 

(1) (WXDE) je = q(WXED) je +af3(WX(D+E)) je - 1 5qW +1 (WX(D+E)) J _ 1/ . 

(2) f3(WXD) je = 5{WXE)^ + af3(WX) j/ _ 1 - -ySq^WX)^^. 

Proof. We claim the following: if (1) is true, then for any word Y in Z)'s and £"s, 
(WXDEY)ji is equal to 

q(WXEDY) je + a(3(WX(D + E)Y) jt - -fdq m+lYl+1 (WX(D + E)Y) j _ 1/ . 

To prove the claim, let Y be any word in D and E. Then (WXDEY)j£ is equal to: 

kt 

i.k 

= q(WXED) lk Y ljke + a/3(WX(D + E)) ik Y ijM - 7 <5g |x|+1 ^(WX(D + E))^ hk Y ijkl 

i,k ik 

= q(WXEDY) je + a/3(WX(D + E)Y) je - 7<5g |x|+|y|+1 (W/X(L> + E)Y)^ U . 



We could have defined matrices D a , D$, Ep, Ej so that they have this combinatorial interpre- 
tation, and then set D = D a + D$ and E = Ep + E 1 . 
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To deduce the final equality above, we applied Lemma 16.71 to the last term. 

Now note that if we take the equation of the claim, and sum over all j and £, 
then we get precisely (I) (since multiplication on the right by V has the effect of 
summing over all indices). And if we take (2) and sum over all j and £, we get (II) 
This completes the proof. □ 

By Proposition 16.8} to prove Theorem 13.4} it is enough to prove the following. 

Theorem 6.9. The following identities hold for all non-negative j and i. 

(1) (WXDE) je = q(WXED) je +aP(WX(D+E)) je - 1 5qW +1 (WX(D+E)) j ^ 1/ . 

(2) P{WXD)jt = 6(WXE)j. ltt + aPiWX)^ - -ySq^WX)^^. 

Note that every term in (1) above is the product of f^ x ^ 2 ~'- l ~ l and a polynomial 
in variables a, 7 and q only. This is easy to see from the combinatorial interpretation 
of the terms: for example, (WXDE)ji is the generating function for certain staircase 
tableaux of size \X\ + 2, with j rows indexed by 5 and £ rows indexed by a/7 (which 
implies that \X\ + 2 — j — £ are indexed by 0. Every term in (2) is also the product 
of <P ^\ x \ +2 ~^~ e and a polynomial in variables 0,7 and q only. 

We will prove Theorem 16.91 by proving both identities simultaneously by induction 
on the degree in and 5 in these identities, and also on the length of the word in 
D and E (equivalently, the size of the corresponding tableaux). Therefore it will be 
convenient for us to make the following notation. 

Definition 6.10. We refer to identity (1) in Theorem \6.9\ as l™ m , where n = \X\+2 
and m = \X\ + 2 — j — £. And we will refer to identity (2) above as 2™ m where 
n = \X\ + 1 and m = \X\ + 2 — j — I. For identity (2), if X = XD (respectively, 
X = XE ), where X is some word in D and E, we will refer to the identity as {2D) 
or 2£>™ m (respectively, (2E) or2EJ m ). 

Our plan for proving Theorem 16.91 is to demonstrate the following implications: 

Proposition 6.11. The following hold for all j, m, and n. 

qririforf < j,m' < m) and l]_ 1>m =► 2D£ TO . 
q&iforj' < j,m' < m) and l£ m _ a => 2E£ m . 
2"w(/or f < J, m' < m) and 2E\ m => l£ m . 

Note that after replacing 2E™ m in the third implication by 2™/^, (for j 1 < j, m! < 
m) and 1" m _ 1 (using the second implication), it is clear that Proposition 16.111 enables 
us to prove identities (1) and (2) in Theorem 16.91 by induction on the sum m + j + n. 
Of course we should also check suitable base cases. 

Lemma 6.12. Theorem \6.iA holds when n < 1 (and j and m arbitrary) and when 
j = (with n and m arbitrary) and when m = (with n and j arbitrary). 

Proof. When n = there is nothing to prove. When n = 1, (1) is vacuous and (2) 
can be proved by inspection. There are only four staircase tableaux of size 1 (single 
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boxes filled with a, (3, 7, or 5), and it is immediate that: 0WD Qj % = a(3Wo^ = at/3, 
(3Di jQ = 6E 0t0 = (35, and 7<Wo,o = SE 0) i. For all other values of j and £, all terms of 
equation (2) are 0. 

Now suppose that j = 0. First consider (1). We must show that (WXDE)j£ = 
q(WXED)j£+a(3(WX(D+E))je, and since j = 0, we should interpret each term as a 
generating function for tableaux which contain no <5's. We will prove this by exhibit- 
ing a bijection for the tableaux corresponding to each term. The term (WXDE)jg 
describes tableaux whose lowest diagonal box is a (3 or 7, with an a in the adjacent 
diagonal box. Consider such a tableau T: if the box just above the (3 or 7 is empty, 
then we map T to the tableau obtained by swapping the entries of the lowest two 
diagonal boxes. The weight of the new tableau is q~ l times the weight of T. Other- 
wise, if the box just above the (3 or 7 contains an a, (3 or 7, then the lowest diagonal 
box must be a (3, and if we delete the lowest diagonal box and the columns second 
from the left (which contains an a with empty boxes above, all weighted by u = 1), 
then we obtain a tableau corresponding to the set a(3(WX(D + E)) ]t . See Figure HI 
It is easy to show that the maps we have described are bijections. 

Now consider j = and (2). We must show that (WXD)ji = a(WX)j ^i. This is 
immediate, as the leftmost column of a tableau from the set (WXD)jg must contain 
an a with empty boxes above, all weighted by u = 1. Deleting the column gives rise 
to a tableau from the set corresponding to a(WX)j t e_i. 

The proofs of (1) and (2) when m = are analogous. □ 
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Figure 4. 



We now prove Proposition 16.111 



Proof of Proposition \6.11\ We start by demonstrating the first implication which 
proves 2_D" m , so we have X = XD, n = \X\ + 1, and m = \X\ + 2 — j — I. 
Note that p{WXD) jt = (3{WXDD)^ = k {WXD) ik D ijkl is equal to: 
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(6) = q8(WXED)^ hi + J2 a(3(WX) hk ^(5(D + E) id ^ 1>£ + A^k-M-i) 

i.k 

i,k 

(7) = qd(WXED) j _ he + a(35(WX(D + E)) j _ 1/ + aPiWXD)^ 

->yPqW + \WX(D + E))^ - jdq^WXD)^^ 

(8) = (6(WXDE)j- lit - af35{WX{D + + ^5 2 q^ +1 {W X{D + S)) j _ 2j /) 

+a/35(W / X( J D + E))j-i t t + a^WIDj^H 
- 7 5 2 g^ l+1 (^X(D + E))j-2,e ~ i5q^(WXD) j ^ 1/ -i 
= 5(WXDE)^ lti + aPtWXD)^! - j5q^ (WXD)^^. 

To deduce (jSj), we use identity 2™~) for z < j and m' < m; note that in the previous 
expression, it was enough to sum over i where i < j. To deduce flBJ) from fl5J), we use 
Lemma I6TT1 to replace by qD^xj-x^t m the nrs t term. To deduce (JSj) from (J7j), 
we apply identity l?_ lm to the term q5(W XED)j_i £. 

For the second implication, we have X = XE. Then /3(WXE)j_i^ is equal to: 

(9) = 5 ^(WXE) ik E id „ 1>k/ 

i . k 

(10) = P(WXD) l+hk E hJ _ hk/ - 

i,k 

(11) = —(WXDE)ji - V + £ + 

+7<5<Z 1 * 1 ^(^Xk^x^Z) + E)ij-i,k-i4 + g^j-i,fc-M_i) 

(12) = Pq-\WXDE) ji -aP 2 q- 1 (WX(D + E)) je -aP(WXE) jti -i 

+/3 7 ^ I ^ I (^X( J D + E))j-i£ + 7 ^l +1 (^ j . v _ 1 

(13) = (/3(Wl££>)^ + a/3 2 g~ 1 (W / X(D + E)) jt - /3j5q^(WX(D + E))^) 

-aP 2 q-\WX{D + E)) jt - a^WXE)^ 
+p-y6qW(WX(D + E))j- U + -ySqW+^WXE)^..! 
= f3{WXED) je - aptWXE)^-! + 'y5qW +1 (WXE) j -x*-i- 

As before, in (jUJ), we may sum over i where i < j — 1. To deduce (ITU]) from (|9l) , we 
use identity 2™~,, for i < j — 1 and m' < m. To deduce (iTTj) from (jTOI) . we replace 
Ei t j-i tk ,£ by q~ l Ei+ij^£ in the first term, using Lemma 16.71 To deduce (lT3"j) from 
(fT2l, we apply identity l\ m _ x to the term pq- 1 (WXDE) j/ . 
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To prove the third implication, note that (WXDE)j£ = k (WXD)i k Eij k £ equals: 
(14) | ^2(WXE)i- ltk E ijkt + a J2( WX kk-xE ijke - ^ WX ^~^ E ^ 



(15) = ^(WXEE) J _ l/ + aY,(WX) hk _ 1 (f3(D + E)^ k _ 1/ + qE 



(3 



i,j,k-l,£-lj 



i.k 

(16) = q(3- 1 5(WXEE) j _ l/ + af3(WX(D + E)) je + aq(WXE) j/ _ 1 

-lSq ]X]+1 J2( WX ( D + E ))i-W ~ P^jSq^iWXE)^-! 

i.k 

(17) = (q(WXED) j£ -aq(WXE) j> ^ 1 + qf3- 1 ^6q^ +1 (WXE) j ^ 1 ) 

+a(5{WX{D + E)) jt + a?(^%i - 7 5g |x|+1 ^(^(D + 

z . k 

-p-i-ySqW+^WXE)^-! 
= q(WXED) je + a/3(WX(D + E)) jt - 7 5g |x|+1 J2( WX ( D + E ))j-U- 

i.k 

In the arguments above, we use 2™~, for i < j and m! < m to deduce ( |14|) : note 
that as before, it was enough to sum over just those % satisfying % < j. To go from 
(JT4"|) to ffT5l) . we replace in the first term by qEi_\^_x, k ,t- To go from f|T6|) to 
(fTTD, we use identity 2££ TO . □ 

6.4. Applications. Once we have a solution to the Matrix Ansatz, it is easy to 
express physical quantities in terms of matrix products [18]. Set C = D + E. 

The partition function Z n is written as WC n V, and the average particle number at 
site i, (Tj) n (where the bracket indicates the average over the stationary probability 
distribution) is written as 

, , WC^DC^V 

(18) {n) = ? • 

Similarly the two-point function (r^Tj) n is given by 

, . WC i - 1 DC*-*- 1 DC n -iV 

(19) (nrj) = , 

A n 

and the n-point functions are expressed similarly. The particle current through the 
bond between the neighboring sites from left to right, which is defined by J = (ri(l — 
Tj+i) — q(l — Ti)r i+1 ), is simply given by J = -jF^. This expression is independent of 
i, as expected in the steady state. 

Note that the matrices D and E that we have defined in Section 16.11 actually 
satisfy the Generalized Matrix Ansatz, not the Matrix Ansatz. However, we can 
compare quantities computed via the two different Ansatzes using Lemma 17.11 and 
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in particular equation (120]) . Theorem 13.51 now follows from Theorem 16.41 and the 
expressions above for the current and m-point functions in terms of matrix products. 



7. The proof of our Askey- Wilson moment formula 
Before proving Theorem 14.21 we need to prove the following result. 

Lemma 7.1. Let D,E,W,V be a solution to the Ansatz of Theorem 15.4 and let 
D, E, W, V be a solution to the Ansatz of Theorem \5. 11 Let h denote the ratio 
Then if X is a word in D and E, and X is the corresponding word in D and E, then 

WXV = h~ l WXV ] [ Xi. 

i=0 

In particular, if Z n = W[D + E) n V , then Z n = h~ l Z n Y\d=Q hi- 
proof. Let t denote the type of X, and let n = \X\. We use induction on n. By The- 
orem [572] and Theorem 15. II respectively. WXV and WXV compute (unnormalized) 
steady state probabilities of being in state r. Therefore WXV = c n WXV for some 
constant c n that depends on n but not X. We want to show that c n = h^ 1 YYi=o \- 

Since we have assumed that D, E, W, V satisfy the relations of Theorem 15.2] 
-fWDXV - aWEXV = X n WXV. By induction, we conclude that 

n-l 

-iWDXV - aWEXV = \ n WXVh- 1 JJ \. 

But also 

-fWDXV - aWEXV = c n+ll WDXV - c n+1 aWEXV 

= c n+1 (^WDXV - aWEXV) 
= c n+1 WXV, 

by Theorem 15.11 This shows that c n+ i = h^ 1 YK=o w hich completes the proof. □ 



We now prove Theorem 14.2] using some results of 



Proof of Theorera \4.S\ Let Z n denote the partition function from [48], i.e. Z n = 
W(D + E) n V, where D, E, W, V are a solution to the Ansatz of Theorem 15.11 and 
WV = h (see [Ml (4.18)]). Here h is as in Section H Then by (Ml Section 6.1], 

'z + z- 1 + 2" 



c 



Therefore 



dz 
c 4m z 



w ((z + z- 1 )/2) 



q 



z + z 



-1 



+ 1 



II 
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which implies that 
1-q 



dz 



-w 



E 

k=0 
n 

E 

k=0 



4niz 

n\ f dz 



{(z + z- 1 )/*) 



z + z 



-1 



k J J c Arriz 



n 



«; ((z + z _1 )/2) 



z + z 



-1 



Inverting this, we get 



E(-D'-"C) 

n=0 v 7 



Z/r, 



By Theorem 13.41 we know that Z n is the generating function for all staircase 
tableaux of size n. By Lemma 17.11 

n-l 

(20) Z n = h 1 Z n l[(af3- 1 6q l ). 



Therefore 



n=0 ^ 



i=0 

fe^ l - q 

n 



2 / nr=o M-7^ 



□ 



8. Open problems 



8.1. Symmetries in the ASEP. Recall that the ASEP has "left-right," "arrow- 
reversal," and "particle-hole" symmetries, which imply Observation 12.21 

Problem 8.1. For each symmetry above, prove the corresponding identity in Obser- 
vation HOI by describing an appropriate involution on staircase tableaux. 

Proving the second identity in this manner is easy. Namely, define a map i by 
letting l{T) be the tableau obtained from T by switching /3's and <5's, and switching 
a'a and 7's; clearly if wt(T) = a h f3 i2 Y 3 5 i4 q is u ie , then wt(i(T)) = a 13 ^ < *7 il <J ia g <8 u« . 
This plus Theorem 13.41 proves the second identity. It remains to find an involution u' 
proving the first identity (the remaining involution can be constructed by composing 
i' with l). A natural guess is to define l'{T) by transposing T then switching a's 
and 5's, and /3's and 7's. This works when q = u, but not for q ^ u. 

8.2. Lifting the ASEP to a Markov chain on staircase tableaux. 

Problem 8.2. Define a Markov chain on the set of all staircase tableaux of size n 
which projects to the ASEP in the sense of [15], such that the steady state probability 
of a tableau T is proportional to wt(T). Such an approach would give a completely 
combinatorial proof of Theorem \3.J\ (This was done in [15] for 7 = 5 = 0.) 
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8.3. "Birth certificates" for particles. In the ASEP, a black particle enters from 
either the left (at rate a) or from the right (at rate 5). Similarly, a "hole" (or a white 
particle) enters from either the left (at rate 7) or from the right (at rate 0). One 
could imagine defining a more refined ASEP, in which each particle in the lattice has 
attached to it its "birth certificate," that is, the information of whether it entered 
the lattice from the left or from the right. Such an ASEP would be a Markov chain 
on 4 n states (all words of length n in a, /?, 7 and 5), which projects to the ASEP 
upon mapping the letters a and 5 to a black particle, and the letters (3 and 7 to a 
white particle. One could then hope to prove an analogue of Theorem 13.41 as follows: 

Problem 8.3. Fix a lattice of n sites, and let S be the set of all 4 n words of length 
n on the alphabet {a, /3,7, 8}, which we think of as configurations of four kinds of 
particles - two kinds of black particles, labeled a and 5, and two kinds of white 
particles, labeled 7 and (3. Define a Markov chain on S with the following properties: 

• particles labeled a and 7 always enter the lattice from the left, and particles 
labeled (3 and 5 always enter from the right; 

• the Markov chain projects to the ASEP; 

• the steady state probability of state (n, . . . , r n ) is proportional to the generat- 
ing function for all staircase tableaux whose border is (n, . . . , r n ). 

8.4. A combinatorial proof of the relations of the Ansatz. In Section El we 
gave a combinatorial proof that D, E, V, W satisfy relation (III) of Theorem I5.2[ by 
translating it into a statement about tableaux. However, we have not yet found a 
combinatorial proof that D, E, V, W satisfy (I) and (II). 

Problem 8.4. Give a combinatorial proof of relations (I) and (II) of Theorem 5.2. 
We note that when q = t, or one of a, /3, 7, S is 0, the above problem is easy. 

8.5. Specializing our moment formula for Askey- Wilson polynomials. 

Problem 8.5. Show directly that our moment formula recovers already-known mo- 
ment formulas for specializations or limiting cases of Askey-Wilson polynomials. 

9. Appendix: Staircase, permutation, and alternative tableaux 

Definition 9.1. [441 136] A permutation tableau T is a Young diagram (where rows 
may have length ) whose boxes are filled with 's and 1 's, such that each column 
contains at least one 1, and there is no which has simultaneaously a 1 above it in 
the same column and a 1 to its left in the same row. The length of T is the sum of 
its number of rows and columns. 

Definition 9.2. [50] An alternative tableau T is a Young diagram (where rows and 
columns may have length ) whose boxes are either empty or filled with left arrows ^— 
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or up arrows f ; such that all boxes to the left of a «— are empty and all boxes above 
an f are empty% The length of T is the sum of its number of rows and columns. 

See [35j for more information about alternative tableaux. 

Proposition 9.3. There is a bisection between staircase tableaux of size n which do 
not contain any 7 or 5, and: 

(1) permutation tableaux of length n + 1; 

(2) alternative tableaux of length n. 

Proof. We first give a bijection from permutation tableaux to staircase tableaux. 
Define a restricted of a permutation tableau to be a which has a 1 above it in the 
same column. A restricted is rightmost if it is the rightmost restricted in its row. 
If T is a permutation tableau, we replace with a every rightmost restricted 0, and 
replace with a t every 1 which is the highest 1 in its column but is not in the top row. 
We replace every other entry of T by an empty box, and delete the top row (but we 
remember the length of the top row by possibly inserting empty columns to the right). 
The result is an alternative tableau, see Figure and the map can be easily inverted. 



For the second bijection, fix a staircase tableau of size n. For i from 1 to n, if 
the ith diagonal box contains an ft, then delete this entry and the column above it 
(this a will correspond to a vertical step in the south-east border of the resulting 
alternative tableau). Otherwise if the ith diagonal box contains a 0, delete this entry 
and the row to its left (this /3 will correspond to a horizontal step in the south-east 
border of the resulting tableau). Then replace each a with an and each /3 with a 
<— , and discard all the other entries. See Figure El □ 
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Figure 5. From a staircase tableau, to a permutation tableau and an 
alternative tableau 



References 

[1] O. Angel, The stationary measure of a 2-type totally asymmetric exclusion process, J. Combin. 

Theory Ser. A, 113 (2006), 625-635. 
[2] R. Askey and J. Wilson, Some basic hypergeometric orthogonal polynomials that generalize 

Jacobi polynomials, Mem. Amer. Math. Soc. 54 (1985), no. 319. 



3 Actually the alternative tableaux of [50] were defined as Young diagrams with blue, red and 
empty boxes; we define them here using left and up arrows, instead of blue and red boxes, following 



24 SYLVIE CORTEEL AND LAUREN K. WILLIAMS 

[3] J. Baik, P.L. Ferrari, S. Peche, Limit process of stationary TASEP near the characteristic line, 

Comm. Pure Appl. Math 63 (2010), 1017-1070. 
[4] M. Balazs, T. Seppalainen, Order of current variance and diffusivity in the asymmetric simple 

exclusion process, Ann. Math. 171 (2010), no. 2, 1237-1265. 
[5] G. Ben Arous, I. Corwin, Current fluctuations for TASEP: A proof of the Prahofer-Spohn 

conjecture, to appear in Ann. Prob. 
[6] A. Borodin, P. Ferrari, Anisotropic KPZ growth in 2+1 dimensions: fluctuations and covariance 

structure, J. Stat. Mech. Theory Exp. 2009, no. 2, P02009. 
[7] R. Brak, J. Essam, Asymmetric exclusion model and weighted lattice paths, J. Phys. A, 37 

(2004), 1483-4217. 

[8] R. Brak, S. Corteel, J. Essam, R. Parviainen, A. Rechnitzer, A combinatorial derivation of the 

PASEP stationary state, Elec. J. Combin. 13 (2006), R108. 
[9] R. Bundschuh, Asymmetric exclusion process and extremal statistics of random sequences, 
Phys. Rev. E volume 65 031911, 2002. 
[10] A. Burstein, On some properties of permutation tableaux, Ann. Combin. 11 (2007). 
[11] S. Corteel, Crossings and alignments of permutations, Adv. Appl. Math. 38 (2007), no 2, 
149-163. 

[12] S. Corteel, P. Nadeau, Bijections for permutation tableaux, European J. Combin. 30 (2009), 
no. 1, 295-310. 

[13] S. Corteel, R. Stanley, D. Stanton, L. Williams, Formulae for Askey- Wilson moments and 

enumeration of staircase tableaux, arXiv:1007.5174 
[14] S. Corteel, L. Williams, Tableaux combinatorics for the asymmetric exclusion process, Adv. 

Appl. Math. 39 (2007), 293-310. 
[15] S. Corteel, L. Williams, A Markov chain on permutations which projects to the asymmetric 

exclusion process, Int. Math. Res. Not. (2007), article ID mm055. 
[16] J. de Gier and F. H. Essler, Slowest relaxation mode of the partially asymmetric exclusion 

process with open boundaries, J. Phys. A: Math. Theor. 41 (2008), 485002 (25pp). 
[17] B. Derrida, E. Domany, D. Mukamel, An exact solution of a one dimensional asymmetric 

exclusion model with open boundaries, J. Stat. Phys., 69 (1992), 667-687. 
[18] B. Derrida, M. Evans, V. Hakim, V. Pasquier, Exact solution of a ID asymmetric exclusion 

model using a matrix formulation, J. Phys. A (1993), 1493-1517. 
[19] B. Derrida, J. Lebowitz, E. Speer, Shock profiles for the partially asymmetric simple exclusion 

process, J. Stat. Phys. 89 (1997), 135-167. 
[20] E. Duchi, G. Schaeffer, A combinatorial approach to jumping particles, J. Combin. Theory 

Ser. A 110 (2005), 1-29. 

[21] F. Essler, V. Rittenberg, Representations of the quadratic algebra and partially asymmetric 

diffusion with open boundaries, J. Phys A 29 (1996). 
[22] P. Ferrari, H. Spohn, Scaling limit for the space-time covariance of the stationary totally 

asymmetric simple exclusion process. Commun. Math. Phys. 265:1-44, 2006. 
[23] P. Flajolet, Combinatorial aspects of continued fractions, Discrete Math. 32 (1980), no. 2, 

125-161. 

[24] D. Foata, A combinatorial proof of the Mehler formula, J. Combinatorial Theory Ser. A 24 
(1978), no. 3, 367-376. 

[25] G. Gasper and M. Rahman, Basic Hypergeometric Series, second edition, Cambridge University 
Press, Cambridge, 2004. 

[26] L. Gwa and H. Spohn, Bethe solution for the dynamical scaling exponent of the noisy Burgers 

equation, Phys. Rev. A 46, 844-854 (1992). 
[27] M. Ismail and D. Stanton, More orthogonal polynomials as moments, Mathematical Essays in 

Honor of Gian-Carlo-Rota, Birkhauser 1998, 377-396. 



COMBINATORICS OF THE ASYMMETRIC EXCLUSION PROCESS 



25 



M. Ismail, D. Stanton and X. Viennot, The combinatorics of the q-Hermite polynomials and 
the Askey- Wilson integral, Eur. J. Comb. 8 (1987), 379-392. 

K. Johansson, Shape fluctuations and random matrices, Comm. Math. Phys. 209 (2000), no. 
2, 437-476. 

A. Kasraoui, D. Stanton and J. Zeng, The combinatorics of Al-Salam-Chihara g-Laguerre 
polynomials, to appear in Adv. Appl. Math. (2010), doi:10.1016/j.aam.2010.04.008. 
D. Kim, D. Stanton and J. Zeng, The combinatorics of the Al-Salam-Chihara g-Charlier Poly- 
nomials, Sem. Loth. Comb. 42 B54i (2006). 

R. Koekoek, P. Lesky, and R. Swarttouw, Hypergeometric orthogonal polynomials and their 
^-analogues, with a foreword by T. Koornwinder, Springer Monographs in Mathematics, 
Springer- Verlag, Berlin, 2010. 

J. MacDonald, J. Gibbs, A. Pipkin, Kinetics of biopolymerization on nucleic acid templates, 
Biopolymers, 6 issue 1 (1968). 

A. de Medicis, D. Stanton and D. White, The combinatorics of g-Charlier polynomials, J. 
Comb. Th. A, 69 (1995), 87-114. 

P. Nadeau, The structure of alternative tableaux, arXiv:0908.4050 

A. Postnikov, Total positivity, Grassmannians, and networks, arXiv:math/0609764vl. 



J. Quastel, B. Valko. t 1 ^ 3 superdiffusivity of finite-range asymmetric exclusion processes on Z, 
Commun. Math. Phys., 273:379-394, 2007. 

T. Sasamoto, One-dimensional partially asymmetric simple exclusion process with open bound- 
aries: orthogonal polynomials approach, J. Phys. A 32 (1999), 7109 - 7131. 
M. Schreckenberg, D. Wolf, Traffic and Granular Flow '97 (Singapore: Springer), 1998. 
L. Shapiro, D. Zeilberger, A Markov chain occurring in enzyme kinetics, J. Math. Biology 15 
(1982), 351 -357. 

R. Simion and D. Stanton, Octabasic Laguerre polynomials and permutation statistics, J. 
Comp. Appl. Math. 68 (1996), p. 297-329. 

D. Simon, Construction of a Coordinate Bethe Ansatz for the asymmetric simple exclusion 
process with open boundaries, J. Stat. Mech. (2009) P07017. 

F. Spitzer, Interaction of Markov processes, Adv. Math. 5 1970, 246-290. 

E. Steingrimsson, L. Williams, Permutation tableaux and permutation patterns, Journ. Comb. 
Th. A, 114 (2007), 211-234. 

C.A. Tracy, H. Widom, Integral formulas for the asymmetric simple exclusion process, Com- 
mun. Math. Phys. 279 (2008), 815-844. 

C.A. Tracy, H. Widom, Asymptotics in ASEP with step initial condition, Commu. Math. Phys. 
290 (2009), 129-154. 

C.A. Tracy, H. Widom, Total current flucturations in the asymmetric simple exclusion model, 
J. Math. Phys. 50 (2009), 09524. 

M. Uchiyama, T. Sasamoto, M. Wadati, Asymmetric simple exclusion process with open 
boundaries and Askey- Wilson polynomials, J. Phys. A. 37 (2004), no. 18, 4985-5002. 
X. Viennot, A combinatorial theory for general orthogonal polynomials with extensions and 
applications. Orthogonal polynomials and applications (Bar-le-Duc, 1984), 139-157, Lecture 
Notes in Math., 1171, Springer, Berlin, 1985. 

X. Viennot, slides and video from the talk "Alternative tableaux, permutations, and 
partially asymmetric exclusion process," at the Isaac Newton Institute, April 23, 2008, 
http : //www .newt on. ac .uk/webseminars/pg+ws/2008/csm/ csmw04/0423/viennot/. 



26 



SYLVIE CORTEEL AND LAUREN K. WILLIAMS 



Laboratoire d'Informatique Algorithmique: Fondements et Applications, Centre 
National de la Recherche Scientifique et Universite Paris Diderot, Paris 7, Case 
7014, 75205 Paris Cedex 13 France 

E-mail address: corteel@liafa.jussieu.fr 

Department of Mathematics, University of California, Berkeley, Evans Hall 
Room 913, Berkeley, CA 94720 

E-mail address: williams@math.berkeley.edu 



